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Abstract 



We give a theoretical overview of the heavy quark tomography puzzle posed by recent non- 
photonic single electron data from central Au+Au collisions at y/s = 200 AGeV. We show that 
radiative energy loss mechanisms alone are not able to explain large single electron suppression 
data, as long as realistic parameter values are assumed. We argue that combined collisional and 
^ ' radiative pQCD approach can solve a substantial part of the non-photonic single electron puzzle. 

in 

1 Introduction 

\Q , Quark Gluon Plasma (QGP) is a new form of matter, consisting of interacting quarks, antiquarks and 

gluons. If the QGP can be created in Ultrarelativistic Heavy Ion Collisions (URHIC), then a wide 
variety of probes and observables could be used to diagnose and map out its physical properties. 

Measured quenching patterns of pions and r\ mesons ^ already provided a direct evidence for 
the creation of a strongly interacting Quark Gluon Plasma (sQGP) in central Au+Au collisions at 
\/s = 200 AGeV [2]-|Sl- Further, rare heavy quark jets are considered to be excellent independent 
probes of the sQGP [S], because their high mass (m c ~ 1.2 GeV, ~ 4.75 GeV) changes the 
sensitivity of the energy loss mechanisms in a well defined way [Zj- J2j relative to those of light quark 
and gluon jets [2]- |U- Another advantage of heavy quarks jet quenching is that gluon jet fragmentation 
into heavy mesons can be safely neglected. However, one disadvantage of heavy meson tomography is 
that direct measurements of identified high pj_ D and B mesons are very difficult with current detectors 
and RHIC luminosities [13] . Therefore, the first experimental studies of heavy quark attenuation at 
RHIC have focused on the attenuation of their single (non-photonic) electron decay products [141 115] . 

The first preliminary data [14] 115] surprisingly suggest that single electrons with p± ~ 5 GeV 
may experience elliptic flow and suppression patterns similar to light hadrons. It was measured that 
the suppression of non-photonic electrons, which is expressed in terms of the nuclear modification 
factor R e AA (P±) = dN(AA -► e) / (N%%dN (pp -> e)), reaches a value ~ 0.1 - 0.4 at p± ~ 4 - 8 GeV. 
Significant reduction at high p± single electrons suggests sizeable heavy quark energy loss. 

Motivated by these data, in JH] we applied the theory of heavy quark radiative energy loss 0- 
[12] to predict the quenching pattern of single electrons from the decay of high p± open charm and 
bottom hadrons. We showed that because the heavy quark "dead cone" effect jZ] is large - especially 
for bottom quarks - radiative energy loss predictions for R\ uAu are significantly above 0.5 as long as 
realistic parameter values are used. Therefore, the puzzle raised by the non-photonic single electron 
data is whether these data can be explained by the energy loss mechanisms in QGP? 

This proceeding mainly concentrate on a theoretical overview of the heavy quark energy loss 
puzzle posed by the single electron data. We start the proceedings with a brief overview of heavy 
quark production and radiative heavy quark energy loss mechanisms in QGP. We then study the 
bottom contribution to the single electron spectra, and show that this contribution is significant, and 
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can not be neglected in the computation of single electron suppression. We then compute the single 
electron suppression from the radiative energy loss mechanisms, and show that radiative energy loss 
alone leads to a disagreement with the single electron data, as long as realistic gluon rapidity density is 
taken into account. Finally, we concentrate on the collisional energy loss, and show that the inclusion 
of this additional mechanism may lead to a better agreement with the single electron data. 



2 Single electron suppression from radiative energy loss 

In this section we will compute the single electron suppression which comes from the radiative en- 
ergy loss mechanisms, and show that the prediction significantly underestimate the single electron 
suppression as long as realistic values of gluon rapidity density are used. 




production medium energy loss fragmentation decay 

Figure 1: A simplified scheme showing how are the non-photonic single electrons obtained from QGP. 

To compute the single electron suppression, we start from Fig.^ which shows the simplified scheme 
of how the non-photonic single electrons are obtained from QGP. We see that, in order to compute 
the single electron spectra, we need to know the initial heavy quark distributions from perturbative 
QCD, heavy flavor energy loss, heavy quark fragmentation into heavy hadrons, Hq, and Hq decay 
into leptons. The cross section is schematically written as (see I17|): 



Ed 3 a{e 
dp 3 



Ejd 3 a(Q) 
dp 3 



P(Ei -> Ef) ® D(Q -> Hq) ® f(H Q -> e) 



where ® is a generic convolution. The electron decay spectrum, /(Hq — > e), includes the branching 
ratio to electrons. The change in the initial heavy flavor spectra due to energy loss is denoted P(Ei — i 
E f ). 



2.1 Initial heavy quark p± distributions 

One of the main advantages of heavy quarks is their large mass (i.e. Mq 3> Aqcd), which, in principle, 
makes perturbative calculations of heavy quark production possible. 

An extensive study of perturbative heavy quark p±_ distributions can be found in the following 
papers |19[ 120] and references therein. By using these papers we can perturbatively compute and 
compare the charm and bottom p± distributions. To compute the initial heavy quark p±_ distributions 
in central rapidity region (\y\ < 0.5) we used the MNR code |20| . We assume the same mass and 
factorization scales as in Ref. |21j . that is we use M c = 1.2 GeV (M& = 4.75 GeV) for charm (bottom) 
mass. For simplicity, we have concentrated only on bare quark distributions (< k 2 ^_ >= GeV 2 ), and 
the runs were performed by using CTEQ5M parton distributions. Initial p± distributions used in our 
computations are shown on the left panel of Fig. |2j From the left panel of Fig. |2] we see that at low 
momentum region, bottom contribution is negligible compared to charm contribution. On the other 
hand, at higher momentum region these two contributions become approximately the same. This is a 
first indication that bottom contribution may become important in the single electron spectrum. 
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2.2 Radiative heavy quark energy loss 



There are three medium effects that control heavy quark radiative energy loss. These effects are 1) Ter- 
Mikayelian, or massive gluon effect |9l I22j. 2) Transition radiation |23l I24j which comes from the fact 
that medium has finite size and 3) Medium induced radiative energy loss 8 , which corresponds to the 
additional gluon radiation induced by the interaction of the jet with the medium. In ^U] we showed 
that first two effects are not important for the heavy quark suppression, since they lead to a change of 
±0.1 in the charm and bottom Raa- We therefore neglect these two effects in this proceedings, and 
concentrate only on the medium induced gluon radiation spectrum, given by: 



dN ind °Fas L f°° 2qV<iq 2 f dk 2 0(2x(l - x)p 



I 



dx 7T X g J (ifS)2 + ( q 2 + M 2 X 2 + m 2J2 J ((| k | _ | q | )2 + ^2)3/2 (( | k | + | q | )2 + ^2)3/2 

f 9 ,9 9 k 2 - M 2 x 2 -m 2 ) 

Here, k is the transverse momentum of the radiated gluon and q is the momentum transfer to the 
jet. M is heavy quark mass, fi ~ 0.5 GeV is Debye mass, X g ~ 1 fm is the mean free path, L ~ 5 fm 
is assumed thickness of the medium, m g:P = fx/y/2 (Hj is gluon mass in the medium, m g>v ~ Aqcd is 



gluon mass in the vacuum and E = \J p\ + M 2 is initial heavy quark energy. We assume constant 
as = 0.3 in this study. 



2.3 Heavy quark distributions before and after quenching 

By knowing the initial heavy quark p± distributions and the heavy quark radiation spectrum, we are 
able to compute the heavy flavor p± distributions after quenching. For this purpose, we generalized 
the multigluon fluctuation approach from |Hj to the case of finite mass. We here give only the final 
results, while for more details the reader should refer to |101 1161 IT7| . 




2 4 6 8 10 12 14 16 2 4 6 8 10 12 14 16 
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Figure 2: Heavy quark p± distributions before and after quenching are shown on the left and right 
panels respectively. For the right panel, assumed gluon rapidity density is dN g /dy = 1000 [51 US]. 
Dark (light) gray curves correspond to charm (bottom) quarks. 

Figure [21 compare p± distributions for heavy quarks before and after quenching. We see that, while 
before quenching (p+p collisions) charm p± distribution is always larger than bottom pj_ distribution, 
after quenching {Au + Au collisions) bottom pj_ distribution starts to dominate the spectra after 
~ 9 GeV. This is expected, having in mind that bottom loose significantly less energy than charm 
quark (compare dot-dashed curves in Fig- EJ 5 an d it is therefore less suppressed than charm quark. 
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2.4 Radiative energy loss prediction for single electron suppression 
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Figure 3: Single electron pj_ distributions before and after quenching are shown on the left and right 
panels respectively. For the right panel, assumed gluon rapidity density is dN g /dy = 1000. Dark 
(light) gray curves correspond to charm (bottom) quark contribution to single electrons. Black curves 
show total (charm and bottom) single electron pj_ distributions. 



In this subsection we show the single electron p± distributions before and after quenching, which are 
obtained after fragmentation and decay of heavy quark p± distributions from the previous subsection. 
For more details on how we obtained these p± distributions, the reader should refer to |16| I19j. 

Figure El compare p± distributions for single electrons before and after quenching. We see that, 

in the case when quarks are not quenched, bottom contribution to single electron spectrum becomes 

comparable to charm contribution at p± ~ 5.5 GeV. For dN g /dy = 1000 case, the crossover between 

charm and bottom contribution is reduced to p± ~ 3.5 GeV. Therefore, in QGP, electrons in the 

p± ~ 5 GeV region have to be sensitive to both b and c quark quenching. 

i 
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Figure 4: Radiative energy loss prediction for single electron suppression. Assumed initial gluon 
rapidity density is dN g /dy = 1000. Dashed (dotted) curve shows what would be the single elec- 
tron suppression if single electrons would have only charm (bottom) contribution. Full curve shows 
the non-photonic single electron suppression by taking into account both charm and bottom quark 
contributions. 




We can now divide p± distributions before and after quenching to obtain the single electron sup- 
pression that comes from radiative energy loss computations (see Fig. 0J. We see that, for realistic 
values of gluon rapidity density, radiative energy loss predicts small single electron Raa ~ 0.7 ± 0.1, 
which is in disagreement with RHIC single electron data |14| II 5j . One possible solution to this problem 
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is to enhance the gluon rapidity density to a maximal value, which would still fit the lower bound- 
ary on pion suppression data. In |16j . we showed that for dN g /dy ~ 3500 the non-photonic single 
electrons can be suppressed to Raa ~ 0.5 ±0.1. However, such large values of gluon rapidity density 
would violate the bulk entropy bounds. Similarly, in Ref. [TS], it was found that a similarly excessive 
transport coefficient [23], q e ff ~ 14 GeV 2 /fm, was necessary to approach the measured suppression of 
electrons from charm jet decay. This finding raised the question of what is the cause for the observed 
discrepancy between theoretical predictions and experimental results. 

3 Collisional energy loss as a solution to the problem? 

Recent studies HZ! suggested that one of the basic assumptions that pQCD collisional energy loss 
is negligible compared to radiative [2H] may be incorrect. In |261 I27| it was shown that, for a range 
of parameters relevant for RHIC, radiative and collisional energy losses for heavy quarks were in fact 
comparable to each other, and therefore collisional energy loss can not be neglected in the computation 
of jet quenching. This result came as a surprise because from the earlier estimates |28| - |31| . the typical 
collisional energy loss was erroneously considered to be small compared to the radiative energy loss. 

However, the computations |26|-|31| were done in an infinite QCD medium, while the medium 
created in URHIC has finite size. A recent paper by Peigne et al. |32j is the first study that made an 
attempt to include finite size effects in the collisional energy loss. This work suggested that collisional 
energy loss is large only in an ideal infinite medium case 1 , while finite size effects lead to a significant 
reduction of the collisional energy loss. However, this paper did not completely separate collisional 
and radiative energy loss effects. Consequently, it remained unclear how important are the finite size 
effects on the collisional energy loss. 
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Figure 5: Fractional collisional energy loss is shown as a function of thickness of the medium for charm 
and bottom quark jets (upper and lower set of curves respectively). Full curves correspond to finite 
medium case, while dash-dotted curves correspond to infinite medium case. Initial momentum of the 
jet is 10 GeV. We assume constant as = 0.3. 

Therefore, it became necessary to consistently compute (only) the collisional energy loss in a finite 
size QCD medium, and see whether the collisional energy loss should be taken into account in the 
computation of jet quenching. In [HHj we provided a detailed study of the th order collisional energy 
loss in a finite size QCD medium created in URHIC. Contrary to |32j . we find that a finite size medium 
does not have a large effect on the collisional energy loss, as shown in Fig. 

In the case of an infinite QCD medium, the collisional energy loss per unit length dE e i jdL is computed by assuming 
that the jet is produced at xq — oo. The energy loss for a finite size medium is than (simplistically) calculated by 
multiplying this dE e i jdl with the thickness L of the medium. 
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Figure 6: Comparison between collisional and medium induced radiative fractional energy loss is shown 
as a function of momentum for charm and bottom quark jets (left and right panels respectively). Full 
curves show the collisional energy loss, while dot-dashed curves show the net radiative energy loss. 
Assumed thickness of the medium is L = 5 fm and A = 1 fm. 

Comparison between collisional and medium induced radiative energy loss [SSI is shown in Fig. |BJ 
We see that collisional and medium induced radiative energy losses are comparable. Therefore, con- 
sistently with the claims in Refs. |2fi| I27| . we see that collisional energy loss is important, and has to 
be included in the computation of jet quenching. 

The results from Fig. Elare still not included in the computation of jet quenching. However, in 
we computed the single electron suppression by taking into account medium induced radiative energy 
loss together with two different computations (TG [23 and BT [301 ) °f collisional energy loss done in 
an infinite QCD medium (see Fig. [7J). We see that with collisional energy loss, it is possible to reach 
single electron Raa ~ 0.5 ±0.1 even at dN g /dy = 1000. While the preliminary data suggest that 
Raa < 0.4, these predictions are consistent with the data within present experimental and theoretical 
errors. Therefore, we may conclude that combined collisional and radiative pQCD approach may be 
able to solve a substantial part of the non-photonic single electron puzzle. 

We note that collisional energy loss in a finite size QCD medium falls between the two different 
computations j2Hl EHJ used in (for more details see [SB])- Therefore, we expect that the single 
electron suppression results - computed with finite size collisional energy loss [SSI " should be inside 
the middle yellow region presented in Fig. 

We also note that (see ^DEl); an y proposed solution of this puzzle must also be consistent with the 
extensive pion quenching data . In |17j we show that the simultaneous inclusion of path fluctuations 
together with radiative and collisional energy loss makes it possible to satisfy R\a < 0.5 ±0.1 without 
violating the bulk dN g /dy = 1000 entropy constraint [HjllSj and without violating the pion quenching 
constraint R\ A « 0.2 ±0.1 now observed out to 20 GeV. 

Finally, note that in this study we used constant coupling as = 0.3. However, the radiative and 
elastic energy losses depend on the coupling, AE rad tx a s and E el oc a| . Future calculations will have 
to relax the current fixed as approximation and allow it to run. Additionally, the running coupling in 
collisional and radiative processes may be different, since these processes might probe different energy 
scales. 

4 Conclusion 

In this proceeding we applied the theory of heavy quark energy loss to non-photonic single electron 
suppression. We showed that bottom quark contribution can not be neglected in the computation of 



6 








p T (GeV) 

Figure 7: The suppression factor, Raa{P±), of non-photonic electrons from decay of quenched heavy 
quark (c+b) jets is compared to PHENIX J3j and preliminary STAR data ^Hj data in central Au+Au 
reactions at 200 AGeV. Assumed initial gluon rapidity density is dN g /dy = 1000. The upper yellow 
band from ^1 takes into account radiative energy loss only, using a fixed L = 6 fm; the lower yellow 
band includes both collisional and radiative energy losses as well as jet path length fluctuations [T7] . 
The dashed curve shows the electron suppression using radiative and TG j^H] collisional energy loss 
with bottom quark jets neglected. Figure adapted from |17j . 



single electron spectra. Additionally, we showed that the recent single electron data lead to significant 
discrepancies with theoretical predictions based only on radiative energy loss, as long as realistic 
values of gluon rapidity density are taken into account. Finally, we introduced the collisional energy 
loss mechanisms, and showed that combined collisional and radiative pQCD approach may be able to 
solve a substantial part of the non-photonic single electron puzzle. 
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